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Abstract 

We construct a four-generation MSSM with rank-4 Yukawa matrices from intersecting D6 
branes on a T 6 / (Z 2 x Z 2 ) orientifold. The Yukawa matrices obtained provide an example of 
Flavor Democracy (FD), where the Yukawa couplings are all nearly equal. Mass hierarchies may 
then be generated by slight perturbations away from FD. We find that it is possible to obtain 
hierarchical masses for the quarks and leptons of each generation and mixings between them. In 
addition, the tree-level gauge couplings are unified at the string scale. Finally, we also construct 
similar models with one, two, and three generations in which the rank of the Yukawa matrices is 
equal to the number of generations in each model. 



1. Introduction. The main challenge of string phenomenology is to exhibit at least one 
string vacuum that describes the physics of our universe in every detail. Despite progress in this 
direction, thus far this goal remains far from achieved. In the past decade, a promising approach 
to model building has emerged involving compactifications with D branes on orientifolds (for 
reviews, see [TJ |2j [31 H]). In such models chiral fermions — an intrinsic feature of the Standard 
Model (SM) — arise from configurations with D branes located at transversal orbifold/conifold 
singularities [5] and strings stretching between D branes intersecting at angles [6j [7] (or, in its T- 
dual picture, with magnetized D branes [Hi El [TO]). A number of non-supersymmetric intersecting 
D-brane models have been constructed that strongly resemble the SM. 

However, non-supersymmetric low-energy limits of string theory suffer from internal incon- 
sistencies of noncanceled NS-NS tadpoles, yielding models that destabilize the hierarchy of 
scales [llj . A resolution of these issues necessarily requires M = 1 supersymmetry. The first 
semirealistic models that preserve the latter were built in Type IIA theory on a T 6 / (Z2 XZ2) orien- 
tifold [HI [13]. Subsequently, intersecting D-brane models based on SM-like, Pati-Salam [14], and 
SU(5) [15] gauge groups were constructed within the same framework and systematically studied 
in Refs. [161 El HI HH] • The statistics of 3- and 4-generation models was studied in [201 HI] ■ Phe- 
nomenologically interesting models have also been constructed on a T 6 /(Z 6 ) orientifold [22] ■ In 
addition, several different models with flipped SU(5) [23] have been suggested within intersecting 
D-brane scenarios [241 125] , as well as models with interesting discrete-group flavor structures [26] . 

Although F-theory model building [27] has received a great deal of attention lately, inter- 
secting D-brane models remain of great interest. While F theory is useful for overcoming the 
problem of phenomenologically necessary but perturbatively forbidden Yukawa couplings, one 
can also evade this problem in intersecting D-brane models by utilizing the Pati-Salam gauge 
group SU(4) x SU(2)i x SU(2)r [IB] , or by incorporating nonperturbative effects such as D-brane 
instantons [28 ] [29] , 130]. Furthermore, the constraints that must be satisfied to construct globally 
consistent intersecting D-brane models are well understood (see, e.g., [3], [3TJ ) . Thus, intersect- 
ing D-brane constructions offer an exciting avenue for model-building, particularly for building 
SM-like and left-right symmetric models based on the Pati-Salam gauge group. 

The SM exhibits an intricate pattern of mass hierarchies and mixings between the different 
generations. One challenge of any string construction is to explain this structure. Within the 
framework of D-brane modeling it was demonstrated that the Yukawa matrices Y a bc ~ exp(— A a b c ) 
arise from worldsheet areas A a b c spanning D branes (labeled by a, b, c) supporting fermions and 
Higgses at their intersections [TJ [32] . This pattern naturally encodes the hierarchy of Yukawa 
couplings. However, for most string constructions, Yukawa matrices are of rank one. In the case 
of D-brane models built on toroidal orientifolds, this result can be traced to the fact that not 
all of the intersections at which the SM fermions are localized occur on the same torus. To date 
only one three-generation model is known in which this problem has been overcome [33j , and for 
which one can obtain mass matrices for quarks and leptons that nearly reproduce experimental 
values. Additionally, this model exhibits automatic gauge coupling unification at the string scale, 
and all extra matter can be decoupled. It should be commented that the rank-1 problem for 
toroidal models can also potentially be solved by D-brane instantons [281 EH [35] . However, the 
conditions for including these nonperturbative effects are very constraining, and at present there 
are no concrete realizations in the literature in which all constraints may be satisfied. 

Although present high-energy experimental data supports just three generations of chiral 
fermions, a fourth generation remains viable as long as the mass of the extra neutrino v' is 
larger than \Mz, and the fourth-generation charged-fermion masses m t ', my, and m T i lie in the 
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correct mass ranges to avoid constraints from direct searches and precision electroweak mea- 
surements [36]. Recent WMAP7 [37] analysis points towards a higher (than three) number of 
relativistic neutrino species, N eS = 4.34+Q;gg, and explanations have been proposed based on 
sterile neutrinos with masses at sub-eV scales [38] (see also Refs. [39] for earlier discussion of 
astrophysical implications of the fourth generation). Furthermore, the existence of a fourth gen- 
eration can introduce additional CP- violating effects [401 SD H2J US] and can have a strong effect 
on Higgs boson searches at colliders |44j . While not directly related to our present focus, these 
considerations suggest possible manifestations of a fourth generation of fermions. 

Interestingly, the principle of Flavor Democracy (FD) [15], where the Yukawa couplings for 
quarks and leptons are all nearly equal, appears to favor the existence of four generations of 
chiral fermions. In this scenario, the observed mass hierarchies are generated as a by-product 
of slight variations away from FD, which may result from the internal geometry of the string 
construction. For example, the down-type quark masses can be generated naturally in the FD 
approach by a Yukawa matrix which is nearly rank 1, while up- type quark masses can be generated 
by a Yukawa matrix which is nearly rank 2, corresponding to near degeneracies in the Yukawa 
couplings. Furthermore, in the FD approach a seesaw mechanism is not necessary to obtain 
small neutrino masses, and the scenario of three light neutrinos with a heavy fourth is naturally 
obtained by considering small perturbations away from FD such that the Yukawa matrix for 
neutrinos is nearly rank 1. This construction is desirable in intersecting D-brane models since 
a Majorana mass term for right-handed neutrinos is perturbatively forbidden and can only be 
obtained nonperturbatively with D-brane instantons. In fact, all Yukawa couplings for quarks 
and leptons can be allowed perturbatively in intersecting D-brane models with a Pati-Salam 
gauge group, and in such models a standard Majorana term of the form Wm — MNN may 
not even be generated through D-brane instantons since it is forbidden by U(1)b-l, which is 
gauged at the string scale. If light, sterile neutrinos turn out to exist, it may be an indication 
that nature does not utilize the seesaw mechanism to achieve tiny neutrino masses. Finally, FD 
may reduce the amount of fine-tuning necessary to obtain hierarchical Yukawa couplings in these 
models since it is natural for the Yukawa couplings to be nearly degenerate, with the observed 
mass hierarchies being generated by slight departures from this degeneracy. 

In the following, we develop a four-generation model constructed from intersecting D6 branes 
on a T 6 / (Z2 x Z2) orientifold in which one can obtain realistic mass matrices. In particular, 
the Yukawa matrices are rank 4, and one can obtain nontrivial masses for each generation, with 
the mass spectrum over the generations being naturally hierarchical. Additionally, the MSSM 
gauge couplings are unified at the string scale. Also, the hidden-sector gauge groups introduced 
to satisfy the Ramond-Ramond (RR) tadpole cancellation conditions become confining at high 
energies. Finally, we construct similar models with one, two, and three generations of matter in 
which the rank of the Yukawa matrices equals the number of generations. 

2. Rank-Four MSSM. The configuration of D branes must obey a number of conditions in 
order to be a consistent model of particle physics. First, the RR tadpoles vanish via the Gauss' 
law cancellation condition for the sum of D-brane and cross-cap RR-charges [3l [4~6] : 





written in terms of the three-cycles Tr a = (n", /") x (n^ , V% ) x (^3 5 2 ^3) that wrap (nj, mj) times 
the fundamental cycles ([%], [bj]) of the factorizable six-torus T 6 = ]^ =1 T?n. Here, the first two 
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Table 1: General spectrum for intersecting D6 branes at generic angles, where I aa > = 
-^UtMli): and/ a06 = 2 3 -' 3 (-^/2 / 3 + z i n 2 n 3 + n i Z 2 n 3 +n i n 2 Z 3 ) _ M is the multiplicity, and 
as and a a denote the symmetric and antisymmetric representations of XJ(N a /2), respectively. 



Sector 


Representation 


aa 


U(iV a /2) vector multiplet and 3 adjoint chiral multiplets 


ab + ba 


M(f, f ) = I ab = 2-e nli(n^ - nlli) 


ab' + b'a 


M(f, f ) = I* = -2-P nlx«li + nlli) 


aa' + a' a 


M(a s ) = \{ha< - \lao&) ; M(a A ) = \{I aa , + \l a0 &) 



two-tori are rectangular: = (j = 1,2), while the third two-torus can be rectangular (/3 = 0), 
or tilted such that /" = 2m" + n" and (3 = 1. In the T-dual picture the tilt of the third cycle 
[a' 3 ] = [a 3 ] + \\b?\ corresponds to turning on a non-zero NS-NS two-form B field. However, it 
becomes nondynamical under the requirement of its invariance under the orientifold projection 
Q7Z [47J. As a consequence, its flux can admit only two discrete values, resulting in two discrete 
values for (3. Each two-torus possesses the complex structure modulus Xj — I built 
from its radii R^f 1 and R%\ M = 1 supersymmetry, which is favored for reasons of underlying 
consistent low-energy theories of particle physics as well as for stability of D-brane configurations, 
is preserved by choosing the angles between the D-brane stacks and orientifold planes to obey 
the condition [13] 

0i + 02 + 03 = mod 2tt , (2) 

with 0? = arctan(2~ /3j x :/ /°/n") and f3\ t 2 = and (3% = (3. This condition can be written in terms 
of wrapping numbers satisfying the two equations 

x A A a + x B B a + x c C a + x D D a = 0, 
A a /x A + B a /x B + C a /x c + D a /x D < 0, (3) 

where 

A a = —fylll, B a = l\n 2 a n\, C a = n^/^n 3 , D a = n a nll a , 

A a = -nln 2 a nl B a = n l f a ll C a = l l a n 2 a ll D a = l l f a nl (4) 

and x A , %b, %c, an d xd are the complex structure parameters [TT| , where xa = A, xb = 
\ . 2#2+As IxiXii xc = X ■ 2 l3l+l3 ' i IxxXzi x d = A • 2^ +/32 /xix 2 , and A is a positive parameter that 
puts the parameters A, B, C, and D on equal footing. Furthermore, the consistency of the model 
is further ensured by the K-theory conditions [3T| US], which imply the cancellation of the Z2 
charges carried by D branes in orientifold compactifications in addition to the vanishing of the 
total homological charge exhibited by Eq. ([T]). In the present case, nonvanishing torsion charges 
are avoided by considering stacks with an even number of D branes, i.e., N a e 2Z. 

Imposing these constraints, we present the D6-brane configurations, intersection numbers, 
and complex structure parameters of the model in Table [21 and the resulting spectrum in Ta- 
ble [3j with formulas for calculating the multiplicity of states in bifundamental, symmetric, and 
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antisymmetric states shown in Table [U Models with different numbers of generations may be 
obtained for different values of the wrapping number n g as well as the third-torus tilt parameter 
0. The observable sector of the models then has the gauge symmetry and matter content of an 
(N g = 2 1_ ^n 3 )-generation SM with an extended Higgs sector. The extra matter in the models 
consists of matter charged under the hidden-sector gauge groups, and vectorlike matter between 
pairs of branes that do not intersect, as well as the chiral adjoints associated with each stack of 
branes. In addition, one has matter in the symmetric triplet representation of SU(2)l as well as 
additional singlets. In order to have just the MSSM at low energies, the gauge couplings must 
unify at some energy scale, and all extra matter besides the MSSM states must become massive 
at high-energy scales. Furthermore, one requires just one pair of Higgs doublets. 

The resulting models have gauge symmetry [U(4)c x U(2) L x U(2)i ? ] obS ervabie x [USp(2 2 ~ /3 (4 — 
n g )) 2 x USp(2 2_/3 ) 2 ]hidden- The hidden sector, as well as the set of complex structure parameters 
required to preserve M = 1 supersymmetry, is different in each of the models with different 
numbers of generations. In particular, in the tilted case two of the hidden-sector gauge groups fall 
out in going from three-generation to four-generation models. The non-Abelian chiral anomalies 
vanish as a consequence of the RR tadpole condition (pQ). The chiral anomalies from the three 
global U(l)s of U(4)c, U(2) L , and U(2) R inducing couplings of the form A a A F 2 , with A and F 
referring to Abelian and non-Abelian gauge fields, respectively, read [49J 

^chiral = 1 J2 N a (Ia(S ~ I a ^)A a A Fj , (5) 

However, these anomalies cancel against the couplings induced by RR fields via the Green- 
Schwarz mechanism [4"9] : 

„4 RR = 8n g A a A (F 2 - F 2 ) + An g A b A F 2 - An g A c A F 2 , (6) 

such that A chnal + A RR = 0. The gauge fields A a of these U(l)s receive masses via linear 
Y2i c^B^ A A a couplings in the ten-dimensional action, with the massless modes given by ker(c"). 
The latter is trivial in the present model, which means that the effective gauge symmetry of the 
observable sector is SU(4) C x SU(2) L x SU(2) fl . 

In order to break the gauge symmetry of the observable sector down to the SM, we split the 
a stack of D6 branes on the first two-torus into stacks a± and a-i with N ai = 6 and N a2 = 2 D6 
branes, and similarly split the c stack of D6 branes into stacks C\ and c<i such that N C1 = 2 and 
N C2 = 2. The process of brane-splitting corresponds to giving a vacuum expectation value (VEV) 
to the chiral adjoint fields associated with each stack, which are open-string moduli. The gauge 
symmetry subsequently breaks down to SU(3)c x SU(2)l x U(1)/ 3B x XJ(1)b-l, where the U(l)/ 3fl 
and U(l) b—l gauge bosons remain massless. The U(1)/ 3JJ x U(l) b-l gauge symmetry may then be 
broken to U(l)y = f U(1) B _ L +U(1) /3J? by giving VEVs to the vectorlike particles with the quantum 
numbers (1, 1, 1/2, -1) and (1, 1, -1/2, 1) under the SU(3) C x SU(2) L x U(l)/ 3/J x U(1) B -l 
gauge symmetry arising from aid x intersections. The full gauge symmetry of the models is then 
SU(3) C x SU(2) L x U(l) y x [USp[2 2 "' 3 (4 - n g )f x USp(2 2 " /3 ) 2 ], with the hypercharge given by 

Qy = \ (Qai — 3<5a 2 — 3Q Cl + 3Q C2 ) , (7) 

where the a-stack charges provide Qb-l and the c-stack charges provide Qzr- 
The gauge coupling constant associated with a stack a is given by 

tol = l«e(/«.)l, (8) 
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Table 2: D6-brane configurations and intersection numbers for a series of Pati-Salam models with 
2 1 ~ f5 n g generations on a Type IIA T 6 / (Z2 x Z2) orientifold, where the tadpole conditions are satis- 
fied without introducing fluxes. The parameter /3 can be zero or one if the third torus is untilted or 
tilted respectively, while the wrapping number n g may take the values 1, 2, 3, or 4. The complete 
gauge symmetry is [U(4) c x U(2) L x U(2) fi ] observable x {USp[2 2 ^(4 - n g )f x USp(2 2 ^) 2 } hidden , 
and the complex structure parameters that preserve M = 1 supersymmetry are xa = %b — 
n g • xc = n g • xd- The parameters (3f give the /3-functions for the hidden-sector gauge groups. 







U(4) c x U(2) L 


xU(2) K xUSp[2 2 ^(4-n 9 )] 2 


x USp(2 2 - 


0)2 














N 


(n 1 ,/ 1 ) x (n 2 ,l 2 ) x (n 3 ,Z 3 ) 


ns 


riA 


b 


b' 


C 


d 


1 


2 


3 


4 


a 


8 


(0,-1) x (1,1) x (1,1) 








2^n g 










1 


-1 








b 


4 


(n 9 ,l) x (1,0) x (1,-1) 


2 1 -*K-i) 


-2 1 -*K - 1) 
















1 





-n g 


c 


4 


(n fl ,-l) x (0,1) x (1,-1) 


-2^K-1) 


2^K-i) 











-1 





n g 





1 


2 2 - /3 (4-?i 9 ) 


(1,0) x (1,0) x (2*0) 


xa — x b — 


n g ■ x c = n g ■ x D <H> xi z 


= rig, X2 




1, 


X3 -- 


= 2* 


2 


2 2 " /3 (4-n s ) 


(1,0) x (0,-1) x (0,2*) 




ft = 


-3, ft 




-3 












3 




(0,-1) x (1,0) x (0,2*) 




PI 


= -6 + 


n 9 












4 


2 2-0 


(0,-1) x (0,1) x (2*,0) 




Pi 


= -6 + 


"< 















where f a is the holomorphic gauge kinetic function associated with stack a, given [3], [TT] in terms 
of NS-NS fields by: 

f a = [nf n£n£ s - 2^< Z£ w 1 - 2 _ ^n^ Z? Zf u 2 - n« I* 1% u 3 ] , (9) 

where K a — 1 for SU(A^ a ) and K a = 2 for USp(2iV Q ) or SO(2iV a ) gauge groups. The holomorphic 
gauge kinetic function associated with SM hypercharge U(l)y is then given by taking a linear 
combination of the holomorphic kinetic gauge functions from all of the stacks [50] : 

fY = \f ai + \{fa 2 +f Cl +fc 2 ). (10) 

Note that in Eq. ([9]), the four- dimensional dilaton s and complex structure moduli u l refer to the 
supergravity basis. These moduli must be stabilized, and gaugino condensation of the effective 
Veneziano-Yankielowicz Lagrangian [51] provides an example of such a mechanism [52]. Gaugino 
condensation in the hidden sectors can play an important role in moduli stabilization, and it 
might provide a top-down reason why three generations is preferred over four. 

From the complex structure parameters, the complex structures U l are determined to be 

U x = n g -i, U 2 = i, U 3 = -f3 + i. (11) 

The dilaton and complex structure moduli are then given in the supergravity basis bjQ 



'4 



(2/% 9 ) 1/2 2tt \n g 



2* 



1/2 



1 See, e.g., footnote 5 of Ref. [50] for the relation between these and complex structures U l 
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Table 3: The chiral and vectorlike superfields of the model, and their quantum numbers under 
the gauge symmetry U(4) c x U(2) L x U(2) R x USp[2 2 " /3 (4 - n g )} 2 x USp(2 2 " /3 ) 2 . 





Multiplicity 


Quantum Number 


Q 4 


Q2L 


Q2R 


Field 


ab 




(4,2,1,1,1,1,1) 


1 


-1 





F L (Q L ,L L ) 


ac 


2^n g 


(4,1,2,1,1,1,1) 


-1 





1 


Fr(Qr, Lr) 


al 


1 


(4,1,1,^,1,1,1) 


1 








X a i 


a2 


1 


(4,l,l,l,iV 2 ,l,l) 


-1 










62 


1 


(l,2,l,l,lV 2 ,l,l) 





1 







64 


n g 


(l,2,l,l,l,l,iV4) 





-1 





X M 


cl 


1 


(1,1,2,^,1,1,1) 








-1 


X C 1 


c3 


% 


(1,1, 2, 1,1, ¥3,1) 








1 


X c3 


&5 


2 1 -' 3 K - 1) 


(1,3,1,1,1,1,1) 





2 





n 




2^k - 1) 


(1,1,1,1,1,1,1) 





-2 





Si 


cs 


2^(n 9 - 1) 


(1,1,3,1,1,1,1) 








-2 


Tr 


CA 


2 1 "/ 3 K - 1) 


(1,1,1,1,1,1,1) 








2 


Sr. 


ab' 


n g 


(4,2,1,1,1,1,1) 


1 


1 









n g 


(4,2,1,1,1,1,1) 


-1 


-1 







ac' 


n g 


(4,1,2,1,1,1,1) 


1 





1 


$1 




n g 


(4,1,2,1,1,1,1) 


-1 





-1 




be 


2n g 


(1,2,2,1,1,1,1) 





1 


-1 


IT I TTl 




2n g 


(1,2,2,1,1,1,1) 





-1 


1 





where 04 = lng s is the four- dimensional dilaton. Inserting these expressions into Eq. (JH1), one 
finds that the gauge couplings are unified as g 2 = g 2 = |#y = g 2 at the string scale Mx, 

with the value of 04 fixed by the value of the gauge couplings where they unify, g 2 (M x ), which 
assumes different values for models with different numbers of generations at Mx = 2.2 x 10 16 GeV: 

g 2 \ Ng=1 (M x ) = 0.275 , g 2 \ Ng=2 (M x ) = 0.358 , g 2 \ Ng=3 (M x ) = 0.511 , g 2 \ Ng=4 (M x ) = 0.895 . 

(14) 

The corresponding string scale is then given by 

Mat = 4?r {-yr ) M Planck , (15) 

where M P i anck is the reduced Planck scale, 2.44 x 10 18 GeV. 

After fixing the value of 4 , one can then determine the values of the gauge couplings for the 
hidden- sector gauge groups at the string scale: 

gW = 2{i - m ™ { 9 Pjl2) ^' ( 16 ) 



where pi = P2 = +1 and pz = Pa = —I- Using the beta-function parameters (3j in Table [2j the 
scale at which each hidden-sector gauge group becomes confining can be calculated: 



A, = M x ■ exp 



2tt 



9 2 (M x )n g ' 



2^r 

(Pi+l)/2 



(17) 



It can then be checked that the hidden-sector gauge groups have sufficiently negative (3j to 
become confining at high-energy scales. To have only one pair of light Higgs doublets, as is 
necessary in the MSSM in order for the gauge couplings to unify, one must fine-tune the mixing 
parameters of the Higgs doublets, specifically by fine-tuning the p term in the superpotential, 
which may be generated via the higher- dimensional operators [33] : 

W D ^SlS> R H k u H l d , (18) 

where yjf are Yukawa couplings, Mst is the string scale, and the singlets S J R are assumed to 
receive string-scale VEVs, while the VEVs of the singlets S l L are TeV-scale. Note that this term 
may only be generated for the models with n g > 1. The exact linear combinations that give the 
two light Higgs eigenstates are correlated with the pattern of Higgs VEVs necessary to obtain 
Yukawa matrices for the quarks and leptons, 

H u , d = V . <d , (19) 



where v l ud = (H l u Thus, at low energies one obtains MSSM-like models with different numbers 
of generations, with gauge-coupling unification ~ 2.2 x 10 16 GeV, and matter charged under the 
hidden-sector gauge groups becomes confined into massive bound states at high-energy scales. 

As has been mentioned, quantities such as gauge and Yukawa couplings depend on the VEVs 
of the closed-string moduli that parametrize the size and shape of the compactified manifold, as 
well as the open-string moduli that parametrize the positions of the D6-branes in the internal 
space, which are associated with the presence of three chiral adjoints in each stack. These VEVs 
should be determined dynamically. While it is not our goal to solve this problem in the present 
work, it should be mentioned that mechanisms do exist by which this can be accomplished. In 
particular, the closed-string moduli can be stabilized in AdS by turning on fluxes in Type IIA |53j . 
In fact, this mechanism has already been demonstrated for the three-generation model [51]. Also, 
gaugino condensation in the hidden sectors can provide another source of closed-string moduli 
stabilization [52] • The open-string moduli may be frozen if the D-branes wrap rigid cycles, a 
possibility that can exist on the T 6 /(Z2 x Z2) orientifold with discrete torsion [55j ESI [57]. An 
example of a four-generation MSSM-like model constructed from D6-branes wrapping rigid cycles 
is given in [19]. We emphasize the possibility of finding a dynamical reason to explain why nature 
chooses a specific number of chiral generations by studying the moduli stabilization problem for 
models with different numbers of generations, such as our mini-landscape of models. 

3. Yukawa Couplings. As one can see from the previous section (note the filler brane stacks 
in Table [3]), only the models with n g < 4 can satisfy the tadpole conditions without introducing 
fluxes. If we take this condition as a constraint, then the only viable models from the top-down 
point of view have N g = 1, 2, 3, 4, 6, and 8. Furthermore, masses may be generated via trilinear 
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couplings for all generations only for those models with a tilted third torus (/3 = 1). If we also 
take this condition as a constraint, then the only viable models are those with N g = 1, 2, 3, and 
4. Additionally, the SU(3)c factor in the SM gauge group is only asymptotically free for SUSY 
models with four generations or less. Thus, the maximum viable number of generations is four. 

The three-generation model has previously been studied in [33]. As mentioned in the Intro- 
duction, this model exhibits rank-3 Yukawa matrices and it is possible to nearly reproduce the 
correct masses and mixings for the three known generations of quarks and leptons. However, 
since a fourth generation has not yet been definitely ruled out, it is worth considering such models 
to see if they can reproduce the observed masses and mixings for the known quarks and leptons, 
while simultaneously satisfying experimental constraints on the fourth generation. In particular, 
the model with n g = 4 and a tilted third torus (/3 = 1) has rank-4 Yukawa matrices, since all 
intersections occur on the first torus, similar to the three-generation model. In contrast to the 
three-generation model, the Yukawa matrices in the four-generation model may potentially be 
flavor democratic, and thus in some sense more natural. Therefore, studying the Yukawa textures 
generated by this model is of great interest. 

A complete form for the Yukawa couplings y{j for D6-branes wrapping on a full compact space 
T 2 x T 2 x T 2 can be expressed as 



Ys 



{ijk} 



h qu Cr a bc $ 



r=l 



6 (r) 



(20) 



where 







S (r) 



K 



(r) 



) = E 



with r = 1, 2, 3 denoting the three two-tori. The input parameters are given by 



(21) 



Sir) 



(r) 



-(r) 



jir) 
L ab 



+ 



jir) k (r) d (r)( /( 



i r )Jr) T ir)Jr) 
T Jo 



rir) 



+ 



jir) 
L bc 



+ 



ab 



he 0* + 



i r )nir) , Ar) Mr) 
ca V h + hb U c 



Jir) 



a' 



d« 

I ri r ) jir) Ar) i 
\ 1 ab 1 bc ±ca I 



ir)irh 
be ta ) 



Ar) Ar) Ar) 
1 ab 1 bc ±ca 



+ 



diry 



(22) 



where the indices i^ r \ j^ r \ and label the intersections on the r th torus, = gcd(I^ } 
and the integer is a function of i^ r \ j^ r \ and corresponding to different ways of counting 



triplets of intersections. The shift parameters et\ ef\ and eh,' correspond to the relative posi- 



,(r) 



tions of stacks a, 6, and c, while the parameters 6a \ \ and 6 C ' ' are Wilson lines associated 
with these stacks. For simplicity, we set the Wilson lines to zero. The brane shifts and Wilson 
line together comprise the open-string moduli, which must be stabilized in a complete model, as 
mentioned previously. For the present work, we treat them as free parameters, as our primary 
interest is simply to see if one can obtain realistic mass matrices; we discussed possible mecha- 
nisms for moduli stabilization in the previous section. Note that although the above formulas 
for the Yukawa couplings are for T 6 = T 2 x T 2 x T 2 , they may be extended to the present case 
T 6 / (Z 2 x Z 2 ) by including all of the orbifold images in the analysis. However, in the present case 
the cycles wrapped by the orbifold images of a stack of D-branes a are homologically identical 



3 (0 
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to the original cycle wrapped by the stack a. In addition, the intersection numbers between the 
cycles defined on the orbifold turn out to be the same as the intersection numbers between those 
on the ambient torus. Thus, the above formulas for T 6 = T 2 x T 2 x T 2 may be used without 
change on T 6 /(Z 2 x Z 2 ). 

We focus only on the first torus, as the Yukawa couplings from the second and third tori only 
produce an overall constant. We label the left-handed fields, right-handed fields, and Higgs fields 
with the indices i, j, and k respectively, which may assume the values 

% G {0, 1, 2, 3} , j G {0, 1, 2, 3} , k G {0, 1, 2, 3, 4, 5, 6, 7} . (23) 

A trilinear Yukawa coupling occurs for a given set of indices that satisfy the selection rule 

i + j + k = mod 4. (24) 

The resulting rank-4 Yukawa matrices, for each of the U, D, and E sectors, have the form 



i=o,4 



^ YoO,j ■ Vj ^01,3+j • ^3+j • V 2 +j 5*03,1+ j ' v l+j \ 

,3+i ' ^3+i Yll,2+j • V 2 +j Yi2,l+j • Vi+j ll3j • Vj 
Y 2 0,2+j • V 2 +j Y 21 

,1+3 ■ V l+j Y 2 2,j ■ Vj *23,3+j • V 3 +j 

\ 5^30,1+i • vi+j Y 3 ij ■ Vj Y 32 ,3+j ■ vs+j Y 3 3 )2 +j • v 2 +j J 



(25) 



where V\. = (Hk+i). These Yukawa matrices are rank 4, although they tend in practice to appear 
numerically close to being matrices of lower rank, due to both the periodic properties of Eq. ( 12TI) 
and the selection rule Eq. ( 124T) . FD Yukawa matrices can therefore arise even when the Higgs 
VEVs are all distinct and of the same order. However, the degree to which FD is broken depends 
upon differences between the Higgs VEVs as well as on the open-string and Kahler moduli, and 
so the detailed values of masses and mixings vary sensitively with the particular values of these 
parameters. Nevertheless, the moduli are not fine-tuned per se, since they turn out to be quite 
different for the U and D/E sectors. Note that the Yukawa matrices for the four-generation 
model with n g = 2, /3 = are at most rank 2, for then the intersection numbers on the first torus 
satisfy gcd(I a b, I ac , I ca ) = 2. Thus, at most two generations can obtain distinct masses in that 
particular variant of the model. 

4. Numerical Analysis. While the model described above supports rank-4 Yukawa matri- 
ces, the question of whether one can in fact obtain a phenomenologically suitable hierarchy of 
masses and mixings (which are somewhat different in the U, D, and E sectors) requires a detailed 
numerical study. In particular, it is of interest to see if the FD approach may be implemented. 
Strictly speaking, the observed masses and CKM elements must undergo RGE evolution to the 
unification scale Mgut i n order to be compared to model predictions. This requirement intro- 
duces two significant difficulties: First, the masses of the fourth-generation fermions, and their 
mixings with the other three, are of course unknown; for purpose of discussion, we simply take 
m t i —my = m T i =400 GeV, a choice that satisfies all current bounds (58]. Indeed, the combination 
of direct observation bounds, electroweak precision tests, and perturbative unitarity constraints 
from heavy-fermion scattering amplitudes place strong constraints on the possible masses, both 
lower and upper limits |59j. Second, a well-known problem of four-generation models, both su- 
persymmetric and not, is the presence of several large Yukawa couplings (due to several fermions 
with electroweak-scale masses) that generate runaway Yukawa couplings at the TeV scale and 
above. While attempts have been made to stabilize the numerical evolution of the RGEs in 
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4-generation models up to Mgut by including new matter fields (e.g., [60]), questions of the 
robustness of such models remain. For our purposes, we simply assume that some supersym- 
metric model exists (see [SI] for a related discussion) in which the Yukawa couplings all remain 
suitably small (from a perturbative viewpoint) to justify RGE evolution up to Mgut ; and take 
the values of the fourth-generation Yukawas at that scale to be 400 GeV divided by the mass of 
the corresponding third-generation fermion, while the Yukawas of the first three generations of 
fermions are those from the 3-generation MSSM RGE evolution to Mqut [621 163] . While such a 
fit is admittedly a hodgepodge from a phenomenological point of view, its purpose is merely to 
provide a proof of principle for the possibility of suitably hierarchical rank-4 Yukawas. 

Specifically, since overall multiplicative factors in each Yukawa matrix Yjj, Ye, Ye remain 
undetermined [Eq. (|20l) ]. we fit to the fermion mass ratios mjjff 10 = rrik/mi in each sector and 
the independent hierarchical CKM elements Vud~V^3> V us . To accommodate the hierarchical 
nature of the masses, we choose ln(m£f 10 ' fit /m^ tl0 ' exp ) as the actual quantities to be fit, which we 
allow to vary by a chosen multiplicative factor (see below) in order to define a unit of \ 2 ■ Since 
the CKM elements are dimensionless, their analogous contributions to \ 2 are l n (l^f t Ai7 CP |)j so 
that the full \ 2 function reads: 



x 2 



£ [ ln «f °- fit Kf °' exp V in(x m )] 2 + J2 MKMPM H*v,*)] 2 > (26) 

ki ij 

where the x m and x Vj ij indicate that a multiplicative discrepancy by x in any observable amounts 
to a unit of x 2 - We chose 10% discrepancy, for mass ratios with indices in 

Eq. (126]) ranging over ki = tc, cu, tf, bs, sd, bb', Tfi, fie, tt', and x V)i j = 1.05 for the diagonal 
ij — ud, cs and x v ^ = 1.1 for the off-diagonal ij = us entries of the CKM matrix. Since all inputs 
are dimensionless, the fit is sensitive only to the ratios of VEVs Vj /v$ , where j = 1, . . . 7, and 
F refers to the appropriate Higgs VEV (that of Hjj or He). The relevant model parameters 
all appear on the first torus; they are given by k' ew^ and the shift parameters e m [i.e., the 
portions of 5^ in Eq. f[2"2"j) that are independent of i, j, and k], with m — 1,2, 3 denoting the 
U, D, E sectors, respectively: 



I ( e (i) _ e W _ 2e « + 2s {1 A (27) 



The numerical inputs and fit values of two sample fits are given in Table HI A few comments 
on the fits are in order: First, we find that the minimization of \ 2 drives the simulation towards 
the vicinity of a stable minimum of the Kahler parameter k ~ 3 even when starting far away 
from it, e.g., /€ init ~ 6. Second, a similar stability of the minimum is observed for the VEVs 
V W>( D \ Third, we find that the VEVs obey the naturalness requirement of being all of the 
same order, i.e., their ratios satisfy Vj U ' D ^ /vq U ' D ^ = 0(1). Numerical simulations demonstrated 
the existence of a large number of solutions with parameters densely populating regions around 
typical ones shown in Table HJ with very small deviations from the computed values. These 
results fall roughly into two classes, with smaller and larger mixing between the third and fourth 
generation of quarks: Note the Vtv entries in Table HI The first set of rows gives typical VEV 
ratios generated by the minimization of \ 2 ■ I n the second set of rows, one sees that the size of 
X 2 is driven by largely by discrepancies from the large fourth-generation fermion masses (~ 400 
GeV), while the mass ratios of the SM fermions agree to a high level of accuracy. The mixing 
between the third and fourth generation is much smaller in the latter case. The values for extra 
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generations that we calculate from the model are consistent with recent fits to flavor physics data 
of Ref. [61] for m t , = 400-600 GeV, \V VV \ = 0.998 ± 0.006 and \V Vb \ = 0.07 ± 0.08. 

The resulting down-quark Yukawa parameters using the input parameters shown in the two 
columns on the left-hand side of Table |4] are 



Y, 



D 



while the up-quark Yukawa matrix is given by 

/ 5.868 14.371 

10.969 5.495 

5.837 14.415 

\ 10.966 5.498 



/ 7.222 


8.015 


7.236 


8.016 \ 


8.001 


8.501 


8.005 


8.500 


7.236 


8.016 


7.222 


8.015 


8.005 


8.500 


8.001 


8.501 / 



(21 



Yr, 



5.837 
10.966 
5.868 
10.969 



14.415 \ 
5.498 
14.371 
5.495 / 



(29) 



As can be seen, these matrices exhibit near degeneracies such that Yjy is close to rank 2, resulting 
in two heavy quarks (£', t) and two light quarks (c, u). is also nearly rank 2, but (it turns 
out) is closer to rank 1 than Yjj, resulting in two heavy quarks (&', b) and two lighter quarks 
(s, d), with b' one being significantly heavier than the others. As discussed earlier, these near- 
degeneracies arise naturally due to the selection rule Eq. ( 1241) and the periodic properties of 
Eq. ( 121]) . However, it must be reiterated that particular numerical values of parameters marking 
the departure from FD are required to obtain these matrices, as seen from the choices of Higgs 
VEVs, the open-string moduli parameters e a ^ c , and the Kahler modulus on the first torus. A 
complete model would of course provide a mechanism for stabilizing the moduli to these values. 
Alternately, some portion of the departure from FD Yukawa matrices may instead arise from 
D-brane instanton-induced couplings, rather than coming entirely from the Higgs and moduli 
VEVs. 

The numerical analysis shows that the fourth-generation masses are naturally larger than 
those of the third generation. Additionally, we find a natural mass hierarchy between each 
generation. However, the model does not account for nonvanishing CKM elements between 
the second and third generation of fermions without significant deviations in mass ratios. The 
nontrivial values for these elements can easily be introduced by incorporating contributions to 
Yukawa matrices from four-point correlation functions [65]. Finally, we note that the above 
analysis provides a nice example of Flavor Democracy, which makes this analysis more natural 
and less fine-tuned compared to that of the three-generation model [33J. 

5. Conclusion. We have constructed a series of MSSM-like models with different numbers 
of chiral fermion generations from intersecting D6 branes on a T 6 /(Z2 x Z2) orientifold. Each of 
these models satisfies all global consistency conditions, including tadpole cancellation, K-theory 
constraints, and conditions for preserving A/" = 1 supersymmetry. For each of the models, we also 
find that the tree-level gauge couplings are unified at the string scale. In addition, for the models 
constructed with one tilted two-torus, we find that the rank of the Yukawa mass matrices for 
quarks and leptons equals the number of generations. Thus, distinct masses for each generation 
and mixings between them may be generated. 

The Yukawa mass matrices for the three-generation model have previously been studied. For 
this model it was found that it is possible to generate mass hierarchies and mixings that nearly 
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Table 4: Experimental input and sample fit values to the model. Note that we have fixed the 
values of the fourth-generation masses to be 400 GeV. 



«' = 3.2452096, 
e 2 = -0.090518048, 

x 2 = 



1 

2 
3 
4 
5 
6 
7 



k£ 



tc 
cu 
ft 
bs 
sd 
b'b 

Tfl 

ue 
r'r 



ud 
us 
cd 
tb 

tv 

t'V 



ei = 0.203163427 
e 3 = -0.090479810 
111 



(D) I (D) 



1.1484243 
0.9603816 
1.2306571 
1.2967893 
1.1462177 
1.4202164 
1.2432859 



m 



_ ratio, fit 
M 



246.1185 
290.6168 
2.7551 
23.7721 
20.0002 
93.7180 
23.7466 
211.1299 
93.7182 



\V; 



fit I 



0.9763 
0.2164 
0.2164 
0.9815 
0.1916 
0.9815 



(u) , (U) 



1.2146412 
0.7789428 
1.4632390 
0.7500008 
2.6474058 
0.8490424 
2.2181503 



rn 



ratio,cxp 
kt 



247.50 
290.60 
2.2857 
35.700 
19.860 
80.000 
21.830 
211.10 
225.1 



ycxp 
ij 



0.9754 
0.22050 
0.22030 
0.9995 



K> = 



62 = 



3.3338221, 
-0.078384580, 

x 2 = 



ei = 0.210962338 
e 3 = -0.079039780 
142 



1 

2 
3 
4 
5 
6 
7 



k£ 



tc 
cu 
ft 
bs 
sd 
b'b 

Tfl 
\1C 

r'r 



ud 
us 
cd 
tb 
tb' 
t'b' 



(D) I (D) 



1.1638620 
0.9282357 
1.2033923 
1.3064203 
1.1390428 
1.4979862 
1.2831043 



m 



ratio, fit 
kl 



247.1696 
290.5549 

5.4376 
35.6756 

19.8532 
134.1494 
21.8330 
211.0802 
134.1025 



IK 



fit! 



0.9759 
0.2185 
0.2185 
0.9981 
0.0610 
0.9981 



(U) I (u) 
'Jj / v o 



1.1950213 
0.7697829 
1.1610830 
0.7459922 
1.2873315 
0.8369313 
1.2145519 



m 



ratio, cxp 
kl 



247.50 
290.60 
2.2857 
35.700 
19.860 
80.000 
21.830 
211.10 
225.1 



T7-CXP 

ij 



0.9754 
0.22050 
0.22030 
0.9995 
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match those that are observed [33J. In the present work, we studied the Yukawa matrices of the 
four-generation model and found that accomplishing the same for the known three generations 
while simultaneously satisfying constraints on fourth generation fermions may be possible. This 
conclusion, of course, comes with the caveat that the problem of the evolution of large Yukawa 
couplings can be ameliorated. Finally, the obtained Yukawa matrices provide a nice implemen- 
tation of Flavor Democracy, and as such seem somewhat less fine-tuned in comparison to the 
three-generation model. In particular, the up-type Yukawa matrix is almost rank 2, resulting 
in two quarks with large masses and two with small masses. The down-type Yukawa matrix is 
nearly rank 1, resulting in one quark with a large mass, and three quarks with smaller masses. 

As commented earlier, the only model of the type under study for which it is possible to 
generate masses and mixings for each generation are those in which the third two-torus is tilted. 
For this subset of models, the maximum number of generations that can be accommodated while 
simultaneously satisfying the tadpole constraints (without including supergravity fluxes) is four. 
Furthermore, it is known that the maximum number of generations in a supersymmetric model 
for which QCD is asymptotically free is also four. 

Although disfavored, four-generation models are still presently viable phenomenologically. It 
is interesting that one can construct such a "realistic" model with four generations. At present, 
the reason that our universe seems to contain only three generations is unknown. Experimentally, 
the Large Hadron Collider (LHC) should be able to determine the question of the existence of 
the fourth generation definitively. From the top-down point-of-view, string theory does not yet 
appear to uniquely determine the number of generations. However, it may be possible to find 
dynamical reasons for singling out three generations once the moduli stabilization question has 
been completely addressed. We leave this question for future work. 
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by the National Science Foundation under Grant No. PHY-0757394. 
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